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Geoinformation Technologies and Spatial 
Analysis of Carbon Dioxide Transport 
Through Border Line 


Geoinformation technologies and methods of spatial analysis of emissions in the border regions have been 
developed and GIS based software has been created for estimating mass of carbon dioxide (CO) emissions 
that goes through border line. Described mathematical models of processes of CO, emissions in the energy sector 
in the border regions take into account the meteorological data. Spatial analysis of carbon dioxide transport 
processes has been done for Ukrainian — Polish border zone in consideration with wind rose. 


Introduction 


It is very important to have full information about country’s greenhouse gases (GHG) 
emissions for participating in new international obligations concerning reducing and monito- 
ring emissions, verification and fulfillment of economical and ecological commitments. Main 
sources of emissions/absorptions processes are located very irregularly on country territory. 
Analysis of spatial distributions of GHG emissions from different sectors of human activi- 
ties on regional level is useful and valuable for projecting effective nature protecting tools, 
study the ways of reducing emissions and uncertainties of GHG inventories. 

Already existed mathematical models of processes of GHG emissions in the energy 
sector do not consider atmospheric diffusion and influence of wind profile. So development 
of mathematical models that will take into account these factors is relevant task. The aim 
of research is to develop mathematical models of processes of GHG emissions in Energy 
sector in the western regions bordering Ukraine, taking into account gas transportation, to 
be more detailed — models, that will give an opportunity to calculate the amount of carbon 
dioxide that goes through fixed border line. 


Carbon dioxide (CO;) properties 


CO, is a naturally occurring gas in the atmosphere and one of the main greenhouse 
gases as it transmits visible light but absorbs strongly in the infrared and near-infrared. CO, 
is also a pollutant, as it is toxic in higher concentrations: 1 % (10,000 ppm) will make some 
people feel drowsy. Concentrations of 7 % to 10 % cause dizziness, headache, visual and 
hearing dysfunction, and unconsciousness within a few minutes to an hour. 

The behavior of dense gases such as CO) can be quite different than for passive tra- 
cers (with neutral buoyancy) due to density effects. For example, a plume of CO) will spread 
laterally over flat terrain even in the absence of wind because it is denser than the surroun- 
ding air [1]. A concentrated CO. plume will generally hug the ground and cross terrain con- 
tours to reach the lowest elevation. Driven by density effects, a CO2 plume will spread more 
quickly than a neutrally-buoyant gas which spreads only through diffusion, giving the coun- 
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ter-intuitive result that ground level concentrations can drop more quickly for a dense gas than 
for a passive gas tracer under calm conditions [2]. The dispersion of COz is highly depen- 
dent on ambient winds. In weak winds, density (buoyancy) effects dominate the dispersion 
process; in high winds, the gas begins to act more as a passive pollutant, with the ambient 
wind overwhelming density effects. Topographic depressions inhibit the movement of the 
CO? plume under weak winds. At wind speeds of greater than 2 m/s, the plume is quickly 
swept out of the domain and diluted, even for rolling hills of amplitude 50 m [1]. 


Mathematical background 


Mathematical model that consider physical properties of carbon dioxide has been 
built. It involves two cases: windy and calm conditions. 

Data required: 

— meteorological data including annual average wind rose; 

— emission data (mass of carbon dioxide in Gg, source coordinates); 

— border line divided into elementary line type objects and their coordinates. 

At any point and any time moment wind can be represented as a vector having a direc- 
tion and value (speed). Although the wind has three dimensions, it is common only to con- 
sider the horizontal components of the wind. The wind direction is the direction from which 
the wind comes. The wind is considered in meters per second, m/s, and at the vertical height 
of 10 meters about ground. 

The mass of carbon dioxide that goes through border line is directly proportional to 
the emission rate. It is therefore important to use emission rates that are as accurate as pos- 
sible. Calculated emissions in Energy sector from every elementary area object of investi- 
gated regions are used as sources of emission [3]. 

Using digital map the Ukrainian — Polish border line was divided into elementary line 
objects. In order to estimate mass of carbon dioxide that transports through border line there 
were made assumptions [2], [4]: 

— Mass that is emitted from the source is assumed to remain in the atmosphere. None 
of the gas is removed through chemical reaction nor is lost at the ground surface through re- 
action, gravitational settlings, or turbulent impaction. 

— Steady — state conditions: meteorological conditions are assumed to persist unchan- 
ged with time, over the time period of transport from the source to elementary line object 
of border line. 

— The average mass profiles at any distance in the crosswind direction, horizontal (per- 
pendicular to the pass of transport) are well represented by a Gaussian, or normal distribution. 

— Mass of gas that goes through border line is assumed to be direct proportional to con- 
centration on border line. 

— Surface roughness is not considered. 

The Gaussian distribution is used to describe the crosswind and vertical distributions 
that result from turbulent mixing that causes dispersion [4], [5]. The height of this curve is 
described by following function: 

(x-u) 


oat |, (1) 


where is the x position where the center of the distribution occurs. The magnitude of the 


1 
f(x) = exp[— 
o.2n 
peak of the distribution is : . The shape of the distribution in the horizontal, whether nar- 
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row or broad, is determined by the magnitude of standard deviation, o . The principle disad- 
vantage of the Gaussian distribution is that it extends from —oo one side to +oo on the other. 
Real plume spreading will be finite. However, from a practical standpoint, the height of the 
Gaussian distribution is very small beyond the limits of +4o [4]. 

It is assumed that x axis is oriented in the wind direction The y axis is oriented in 
crosswind direction. 

As mass rate of CO, that passes through border line is directly proportional to gas 
concentration, the next formula takes place: 


M(x,y)=K, *c(x,y), (2) 
where M (x, y) is the distribution of mass in crosswind direction, g/m; K,, is the coefficient 


. A 2 5 Z 
which value depends on effective volume of gas plume, m’; c(x,y) is the gas concentration 
, ‘ 5 7 3 
in crosswind direction, g/m’. 
Instantaneous gas concentrations in crosswind direction can be described by formula: 
y? 


2M “a6? 
ayaa aes (3) 
(27) C68: 


where W/,, is the total mass of emissions in g; o, is the standard deviation of concentration 


distribution in the puff in the upwind-downwind direction; o,, is the standard deviation of 


concentration distribution in the crosswind direction, m, at the downwind distance x; o, is 
the standard deviation of concentration distribution in the vertical direction, m, at the down- 
wind distance x [4]. 

The speed of the wind mainly serves to give the download position of the center of 
the puff. Wind speed may influence the dispersion indirectly because the dispersion para- 
meters o,,0,,0, may be functions of wind speed. Much less is known of dispersion in the 


upwind-downwind direction than is known of lateral and vertical dispersion. In general, one 
should expect the o,, to be the same as o,. 
For calculation of o,and o, parameters based on measurements of atmospheric tur- 


bulence over flat plains are used. As it has been observed when air is stable, vertical mixing 
is inhibited, and when there is strong solar heating of the surface, there may be strong con- 
vective activity with large vertical motions. There have been produced empirical results for 
the variation of o,,and o, for six stability classes: A — extremely unstable, B — moderately 


unstable, C— slightly unstable, D — neutral, E — slightly stable, F — moderately stable. 
Guidelines are given for estimating the stability class from wind speed, cloud cover and 
time of day in Table 1. 


Table 1 — Guidelines for determining Pasquill — Gifford stability classes [4] 


: Day with insolation Night 
sia ae St Mod Slioh Overcast > 4/8 <3/8 
speed, m/s rong oderate ight jenrcladdl diana 
2 A A-B B — — 
2-3 A-B B C E F 
3-5 B B-C C D E 
5-6 C C-D D D D 
6 C D D D D 
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Approximate power-law functions (4) can be specified for quasi-instantaneous sour- 
ces corresponding to the Pasquill stability classes [4] (Table 2). 


C= ax’ .o. =x", (4) 
Table 2 — Quasi-instantaneous Power Functions [4] 
Pasquill b Oy Oy d Oz Oz 
stability : 100m | 4km | * 100m | 4km 
A 0,18 0,92 12,45 371 O72 0,76 23,8 393 
B 0,14 0,92 9,69 288 0,53 0,73 15,3 226 
C 0,1 0,92 6,92 206 0,34 0,72 9,4 133 
D 0,06 0,92 4,15 124 0,15 0,70 3,8 50 
E 0,045 0,91 297 85 0,12 0,67 2,6 31 
F 0,03 0,90 1,89 52 0,08 0,64 1,5 16 
“G” 0,02 0,89 1,21 32 0,05 0,61 0,8 8 


If to substitute the valuable c (x, y) (formula (3)) into formula (2): 
a 


2M “262 
E e y . 


(22) CG se: 


M(x,y)=K, * (5) 


The area under curve M(x, y) is assumed to be equal to / ,,, total mass of emissions (see 
Figure 1): 


M, = [M(x,y)dy 


—00 


or 
t 2M 
M,=)|K,* B ’ (6) 
i (22) o,0,0, 
From (6): 
K,=70,0,. 


Final formula for mass distribution in crosswind direction is: 


ee 
M(x,y)=M,, * eee 
20 


Figure 1 — Mass distribution in crosswind direction 
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Mass that goes through elementary line object M,,,(x,,,,) (see Figure 1) can be 


estimated by following formula: 


v2 1 ae 
M gas(%VysV2) = | Me * e ne (7) 


V1 


ae 


ON 
Z 
cA 27 (iy) 
7 " 


4 a \ 


source Xx 


Figure 2 — Mass distribution in case of calm conditions 


In case of calm conditions the average mass profiles at any distance in any direction 
cannot be represented by a Gaussian curve. So it is assumed that mass of gas has even distri- 
bution in all directions. 

Mass of carbon dioxide that passes elementary line object is directly proportional to 
the value of angle a (in radians) that is created by two radii drawn from the ends of line ob- 
ject to source as it is shown in Figure 2. Using cosine rule, here the formula (8) for estima- 


ting rate of mass: 
ate -b? 
arccos| ————— 
2ac 


20 (8) 


M2" (a,b,c)=M , 


gas 


where M el is the mass of gas in calm conditions; a,b are the radii connected source cent- 


re and ends of elementary line object, m; c is the length of elementary line object, m. 


Figure 3 — Moving from coordinate system of digital map to another one with centre 
in the source of emissions and x axis oriented in the wind direction 
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As in practice there is a need to work with big number of sources of carbon dioxide 
emissions and different wind directions, so there were developed formulas below (9) — (10) 
for moving from one coordinate system of digital map to another with center in source of emis- 
sions and x axis oriented in the wind direction. 

Variables in Figure 3: x,y are the coordinates of emission source in “old” coordinate 
system; x’, y’ are the new coordinates with centre in source of emissions and axis x oriented 
in the wind direction; (x;,v,) are the coordinates of ends of line type object in map coordinate 
system; a is the angle between wind direction and x axis; f is the angle between vector, con- 
nected source and receptor, with x axis. 


x, =acos(B—a), y, =asin(B—a@), (9) 


where (x,,,) are the coordinates of ends of line type object in new coordinate system with 


centre in source point and axis x oriented in wind direction. 
The value of angle / can be calculated using formulas (10): 


arctg ead 9X, > XV, > V3 
bee 

7 — arctg oA a <X,V, <3 

a 
2 — arctg oA nae > XV <3 

Xe. 
7m — arctg aa ge <X,V, > Vs 

Pater: 
pry (10) 
9X, = X,Y, > V3 


5 
al ee <y; 
ae Vi Sd; 


0.x, > X,Y, =yV 


2H,X, >X,y, =¥. 


Numerical modeling 


For estimating mass that goes through border line as input data there were used avera- 
ge annual wind rose for western Ukrainian regions (wind speed, and direction with frequen- 
cy), and digital maps with information about carbon dioxide emissions in Energy sector of 
western regions in Ukraine on elementary objects level [1]. 

The tools of geoinformation system (GIS) MapInfo was used. Additional menu was crea- 
ted on the panel, based on developed mathematical model, and programmed on MapBasic. 
It gives opportunities to build border line, select border zone of fixed width with emissions 
sources, and calculate mass that passes through selected border line (see Figure 4). 

For calculations the Ukraine — Polish border line was chose and border zone with 100 km 
width was selected. As a result the digital map with information about mass of carbon dioxi- 
de that goes through border line on elementary object level was created. Figure 5 shows the 
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distribution of mass of carbon dioxide that goes through Ukraine — Polish border line on ele- 
mentary line object level from south to north. 


Carbon dioxide transport 


Figure 4 — MapInfo menu 


Mass of CO?, thousands ton 


Elementary line type objects of border line 


Figure 5 — Mass distribution of COz across border line from south to north 


Conclusions 


Created mathematical model that takes into account meteorological data and is based 
on geoinformational technologies gives the opportunity to estimate rate of mass of carbon 
dioxide that goes through border line. Made estimations for Ukrainian — Polish border line 
and border zone with 100 km width show that only 15 % of CO» emitted in this zone are tran- 
sported through border line from Ukraine to Poland. Creation of common geoinformation 
system for GHG spatial inventory, and formation of georeferenced databases with input da- 
ta and inventory results for Ukrainian — Polish border regions is in perspective. Developed 
methods and geoinformation technologies can be used for estimating mass that transports 
through border line as well. 
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M.FO. JTecue, P.A. byno 

TeonndopmaltvoHuble TexHO.JIOrH WH NpOCTpaHCcTBeHHbIi aHa.IH3 Npoweccos 

TlepenHoca AMOKCHAa yriepoya uepe3 rTpaHHny 

Onncanbt reouHopMallHOHHble TEXHOJIOIMM HM MCTOAbI IPOCTpaHCTBeHHOTO aHaIH3a IMUCCHH WapHMKOBBIX 
Ta30B B IIpHIpaHWyHBIX perMouax Mf CO3{aHoO TIporpaMMuHoe obeciie4eHHe JIA YHCIICHHOTO MOJesIMpOBaHHA TIPOLeCCOB 
llepeHoca WuoKcuya yruiepoya uepe3 rpanuny. MpeqioxKenHbie MaTeMaTH4ecKHe MOJeJIM MpoweccoB IMHCCHU 
yIIeKMCJIOFO Ta3a B 9HEPreTHYeCKOM CeKTOpPe 3alaHbIX perHOHOB YkKpauvHbl JIA BEIYMCIeHuA TepeMelleHuA 
aTMOC(PepHbIX MaCC YYHTHIBAIOT METEOPOJIOFM4eCKHe YCIOBHA, a MMeHHO — po3y BeTpos. IIpoctrpaHcTBeHHBI 
aHaJIM3 IMMCCHM yrIeKHCOrO Ta3a OBI CleaH JIA YKpaMHCKO-MOIbCKOM NOrpaHH4HoH MOJIOCEL. 


M.FO. JTecie, P.A. byno 

Teoindopmaniini TexHoorii Ta MpocropoBuii anai3 Nponecis 

Tepenocy AiokcHAYy Byr1elIo Yepe3 KOPAOH 

OnucaHo reoindopMaliiiHi TeXHOJOrii Ta METOAM MpocTopoBoro aHai3y eMIciit HapHMKOBHX Ta3iB B TIPHKOPAOHHHX 
perioHax Ta CTBOpeHO NporpaMHHit 3aciO [VIA UHCIOBOTO MOJeIIOBAHHA IIpOleciB MepeHocy AIOKCHY ByTsIeLIo 
Yepe3 JIHIKO KOpAOHY. 3anponoHoBaHi MaTeMaTHYHI MOJes MpOLeciB eMicli ByIeKMCIOrO ra3y B eHepreTHYHOMy 
CeKTOpi 3aXiJHUX perionis YkpaiHH JVId OOUMCICHHA TepeMILeHHA AaTMOC(PepHUX Mac BpaXOBYIOTb MeTeOpOJIOT1YHI 
YMOBH, a caMe — po3y BiTpiB. [IpocropoBuii aHali3 eMicili ByrIeKHCOro Tra3y 3poOseHO TIA yKpalHCbKO-MOJIBCBKOI 
IIPHKOPAOHHOi CMyTH. 


Cmamba nocmynua 6 pedaxyuro 21.06.2010. 
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